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ABSTRACT 

Measurements of cosmic rays by experiments such as ATIC, CREAM, and PAMELA in- 
dicate a hardening of the cosm ic -ray energy spectrum at TeV energies. In our recent work 
dThoudam & Horandelll2012al) . we showed that the hardening can be due to the effect of 
nearby supernova remnants. We showed it for the case of proton and helium. In this paper, 
we present an improved and detailed version of our previous work, and extend our study to 
heavier cosmic -ray species such as boron, carbon, oxygen, and iron nuclei. Unlike our previ- 
ous study, the present work involves a detailed calculation of the background cosmic rays and 
follows a consistent treatment of cosmic -ray source parameters between the background and 
the nearby components. Moreover, we also present a detailed comparison of our results on the 
secondary-to-primary ratios, secondary spectra, and the diffuse gamma-ray spectrum with the 
results expected from other existing models, which can be checked by future measurements 
at high energies. 
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1 INTRODUCTION 



Recen t measure ments of cosmic r ays by the ATIC jPanov et al 



20071). CREAM JYoon et alj[20nh . and PAMELA JAdriani et al 



20111) experiments have indicated the presence of hardening in 



the energy spectra of protons, helium and heavier nuclei at TeV 
energies. The observed hardening does not seem to be in good 
agreement with the standard theoretical predictions. According 
to th e standard theory of cosmic-ray acceleration (IBelll 119781 : 
iBlan dford & Eichler 1987|) and the nature o f their transport in the 
Galaxy (see e.g., G inzburg & Ptuskinlll976r) . the cosmic-ray spec- 
trum should follow a single power-law at least up to ~ 3 PeV, the 
so-called "cosmic-ray knee". 

Different explanations for the spectral hardening have been 
proposed. Most of these explanations suggest either hardening in 
the cosmic-ray source spectrum or changes in the propagation 
properties in the Galaxy. Possible scenarios that can produced a 



hardened source spectrum have been sugg este d in|B iermann et al. 
l20ld lOhira et aLllioTTi I Yuan et al.l l201 1 and Ptuskin et al. 2012, 
and models bas ed o n propagation effe cts have been discussed in 
lAve et alj 1 20091 and iTomassettil 1201 2| . Various interpretations of 
the spec tral hardening including the effect of local sources can be 
found in lVladimirov et al.ll2012i 

It is generally believed that nearby sources can affect the prop- 
erties of cosmic rays observed at the Earth, and they might also ac- 
count for the observed spectral hardening. In lErlvkin & Wolfendald 
I2OI2I it was suggested that the hardening can be due to a steep lo- 



E-mail: s.thoudam@astro.ru.nl 



cal component dominating up to ~ 200 GeV/n, and a flatter back- 
ground dominating at higher energies. Another scenario is that the 
high-energy spectrum might be dominating by the local component 
and the low-energy region is dominated by the background. This 
scenari o was proposed in our r ecent wo rk (IThoudam & Horandell 
l2012al . hereafter referred to as IPaper 11) . in which we considered 
nearby supernova remnants within 1 kpc in the Solar vicinity. In 
addition to explaining the spectral hardening, we could also explain 
the observed hardening of the helium spectrum at relatively lower 
energy/nucleon with respect to the protons. Moreover, we showed 
that the hardening may not continue up to few tens of TeV/nucleon. 
In iPaper J . we concentrated only on protons and helium nu- 
clei. Moreover, we did not perform a detailed calculation for the 
background component. The background was obtained by fitting 
the measured low-energy data in the range of (20 — 200) GeV/n. 
Considering that this is the energy region where both the contri- 
bution of the local sources and the effect of the Solar modula- 
tion are expected to be minimum, the background thus obtained 
seems to represent fairly well the averaged cosmic-ray background 
present in the Galaxy. However, the background should be con- 
sistent with the observed data down to lower energies say up to 
su b-GeV o r ^ 1 GeV energies. This consistency was not checked 
in lPaper j . In addition, if we assume a common cosmic-ray source 
throughout the Galaxy, then the source parameters for the back- 
ground component such as the spectral index and the cosmic-ray 
injection power should also be consistent with those of the nearby 
sources. To check whether the spectral hardening can still be ex- 
plained by maintaining such consistencies, a detailed calculation 
should be carried out for both the components together This is the 
main motivation of our present paper A similar study was recently 
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presented in lBemard et al.l2012i The main difference between their 
work and ours is on the treatment of the nearby sources. They as- 
sumed supernova remnants as instantaneous point-like sources in- 
jecting cosmic rays in an energy-independent manner, while we 
consider finite-sized remnants producing cosmic rays of different 
energies at different stages in their hfetime. For distant sources, 
the much simpler energy-independent point-source approximation 
may represent a valid approximation, but for nearby sources, it is 
more realistic to take into account t heir finite sizes and the nature o f 
cosmic-ray injection in the Galaxy ( IThoudam & Horandelll2012b[) . 
In the present work, we will further extend our study to heav- 
ier cosmic-ray species namely boron, carbon, oxygen, and iron. 
We will also present a detailed comparison of our results on the 
secondary-to-primary ratios, secondary spectrum and the diffuse 
gamma-ray spectrum with those expected from other existing mod- 
els. Our model is described in section 2. Then, our calculations for 
the cosmic-ray spectrum from a nearby source and for the back- 
ground spectrum are presented in sections 3 and 4 respectively, and 
in section 5, the various interaction cross-sections and the matter 
density that will be used for our study will be presented. In section 
6, constraints on the contribution of the nearby sources imposed by 
the secondary cosmic-ray spectrum are discussed. In section 7, the 
results on the heavy nuclei are presented, and in section 8, the re- 
sults on protons and helium nuclei are presented. In section 9, we 
present the discussions and compare our results with the predic- 
tions of other models, and we conclude with a short summary in 
section 10. 



2 THE MODEL 

Although there is no direct evidence yet that proves that supernova 
remnants are the major source of Galactic cosmic rays, they repre- 
sent the most probable candidates both from the theoretical and the 
observational point of views. Theoretically, it has been established 
that supernova remnants can produce a power-law spectrum of par- 
tic les up to very high energies with a sp ectral index very close to 
2 JBellll 19781 : fBlandford & Eichleilll987h . This value of the spec- 
tral index agree nicely with the values determined from the obser- 
vations of supernova remnants in radio frequencies (Green 20090 . 
Moreover, if ~ (10 — 20)% of the total supernova kinetic energy 
is channeled into cosmic rays, supernova remnants can easily ac- 
count for the total amount of cosmic-ray energy contained in the 
Galaxy. Observationally, the presence of high-energy particles in- 
side supernova remnants is evident from detections o f non-thermal 
X-rays and TeV 7-rays from a num ber of remnants (IParizot et al.l 
1200a lAharonian et al.ll2006l I2OO8I) . The non- thermal X-rays are 
most likely synchrotron radiations produced by high-energy elec- 
trons in the presence of magnetic fields. The TeV 7-rays are pro- 
duced either by inverse compton interactions of high energy elec- 
trons with ambient low energy photons or by the decay of neutral 
pions, which are produced by the interaction of hadronic cosmic 
rays, mainly protons with the surrounding matter. Irrespective of 
the nature of production, the detection of TeV 7-rays suggests the 
presence of charged particles with energies larger than few TeV in- 
side supernova remnants. 

Based on these observational evidences and our current theo- 
retical understanding, we assume that supernova remnants are the 
main sources of cosmic rays in the Galaxy. We further assume that 
cosmic rays observed at the Earth consist of two components: a 
steady background which dominates most of the spectrum and a 
local component which is produced by nearby sources. The back- 



Table 1. List of supernova remnants with distances < 1 kpc considered in 
our study. 



Name 


Distance (kpc) 


Age (yr) 


Geminga 


0.15 


3.4 X 10^ 


Loopl 


0.17 


2.0 X IQS 


Vela 


0.30 


1.1 X 10" 


Monogem 


0.30 


8.6 X 10" 


G299.2-2.9 


0.50 


5.0 X 103 


Cygnus Loop 


0.54 


1.0 X 10" 


G114.3+0.3 


0.70 


4.1 X 10" 


Vela Junior 


0.75 


3.5 X 103 


S147 


0.80 


4.6 X 103 


HB9 


0.80 


6.6 X 103 


HB21 


0.80 


1.9 X 10" 


SN185 


0.95 


1.8 X 103 



ground is considered to be produced by distant sources which fol- 
low a uniform and continuous distribution in the Galactic disk, and 
the local component is assumed to be contributed by supernova 
remnants with distances within 1 kpc from the Earth. A list of su- 
pernova remnants that will be considered for our study are given in 
Table 1. 

The two cosmic-ray components are treated with different 
propagation models. The background component is treated in the 
framework of a steady state diffusion model, and the local compo- 
nent in a time dependent mod el. For the backg round component, we 
adopt the model described in lThoudamll200a where the cosmic-ray 
diffusion region is taken to be a cylindrical disk with infinite radial 
boundary and finite vertical boundaries ±_ff. For a typical value 
of the radial boundary which is ~ 20 kpc or more, our assump- 
tion of infinite radial boundary represents a good approximation 
for cosmic rays at the Earth. It is because for large radial bound- 
ary, cosmic-ray escape is dominated by escape through the vertical 
boundary and the effect of the radial boundary on the cosmic-ray 
flux becomes negligible. Regarding the size of H, diff erent models 
adopt different values in the range of --^ (2 — 10) kpc dStrong et al.l 
uOlOl) . and for the present study, we assume H = 5 kpc. Further- 
more, both the background sources and the interstellar matter are 
assumed to be distributed in the Galactic plane in an infinitely thin 
disk of radius R. And we take R = 20 kpc for our study. 

For the local component, we assume a diffusion region with 
infinite boundaries in all directions. This assumption is made, con- 
sidering that the spectrum of cosmic rays from nearby sources do 
not depend much either on the radial or the vertical boundaries be- 
cause of their very short propagation time to the Earth relative t o 
the escape times from the Galactic boundaries ( IThoudaml 120071) . 
Strictly speaking, the propagation does depend on the boundary be- 
cause of the dependence of the diffusion coefficient on the vertical 
boundary. But, once the diffusion coefficient is fixed, one can ne- 
glect the dependencies on the diffusion boundaries and determine 
the cosmic-ray flux using infinite boundaries. 



3 COSMIC RAYS FROM A NEARBY SOURCE 

The propagation of cosmic rays from a nearby source is governed 
by the diffusive transport equation, 



dN 
V ■ (DVN) +Q = -Q^ 



(1) 
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The first and the second terms on the left represents the diffusion 
and the source terms respectively, with A'^(r, E, t) representing the 
particle number density, E the kinetic energy/nucleon, r the dis- 
tance from the center of the source and t representing the time. The 
diffusion coefficient D is assumed to be a function of the particle 
rigidity p as D{p) — Dol3{p/po) , where 5 is the diffusion in- 
dex, P — v/c with V denoting the particle velocity and c is the 
velocity of light. For a nuclei carrying charge Ze and mass num- 
ber A, the rigidity can be written as p = APc/Ze, where P de- 
notes the momentum/nucleon of the nuclei. Thus, for the same P 
(~ E for all the energies of our interests here), all heavy nuclei with 
charge Z > 1 diffuse relatively faster than the protons by a factor 
of (A/Z) . In Eq. (1), we neglect the effect of interactions with 
the interstellar matter as the time taken for cosmic rays to reach 
the Earth is expected to be much less than the nuclear interaction 
timescales. Also, we do not include effects which are important 
only at low energies such as the ionization losses, convection by 
the Galactic wind, and the possible re-acceleration by the interstel- 
lar turbulence. With these approximations, we assume that Eq. (1) 
effectively describes the propagation of cosmic rays with energies 
above 1 GeV/n. 

As mentio ned before, an important feature of our model (as 
also in IPaper m is the assumption of finite-sized sources which 
inject cosmic rays in an energy-dependent manner. Our model is 
based on the current basic understandings of diffusive shock accel- 
eration (DSA) theory inside supernova remnants. In DSA theory, 
charged part icles are accelerated each time they cro ss the supernova 
shock front ( lBelllll978l:lBlandford & Eichlerll987h . During the ac- 
celeration process, a major fraction of the particles are advected 
downstream of the shock, while a small fraction escape upstream 
which can be again taken over by the expanding shock. This allows 
repeated crossings of the upstream particles across the shock front. 
Thus, efficient acceleration can be achieved when particles are ef- 
fectively confined near the shock in the upstream region. It is now 
generally accepted that the confinement can be provided by mag- 
netic scattering due to the turbulence generated ahead of the shock. 
Particles can remain confined as long as their upstream diffusion 
length is much less than the shock radius, i.e.. Id "C Ra- The dif- 
fusion length is related to the diffusion coefficient in the upstream 
region Ds and the shock velocity Us as Id = Dg/us- Then, for Ds 
that scales with energy, for instance Ds oc i? in the Bohm diffusion 
limit. Id oc E/us- This shows that high-energy particles can escape 
at early stages during the evolution of the remnant while the shock 
is still strong. For the lower energy particles, they can escape only 
at later stages when the shock becomes weak. With this general un- 
derstanding, but at the same time considering that the actual energy 
dependence of Ds is not clearly understood because of many com- 
plicated processes involved during the acceleration process such as 
magnetic field amplification and back-reaction of the particles on 
the shoc k, the co smic-ray escape time is parameterized in a simple 
form as jPaper M) 



top„ 



1 PV. Expressing Eq. (2) in energy/nucleon, we get 



f-esc^-JtJ — f'sed 



-1/q 



(2) 



i'esc \-^ } — ^i 



V Zepm ) 



-1/q 



(3) 



Eq. (3) shows that for the same E, all heavy nuclei escape at 
relatively early times compared to the protons by a factor of 
{A/Z)^^'°' . This early escape of heavier nuclei in our model was 
the key to explaining the observed spectral hardening of helium at 
lower energy/nucleon with respect to the protons in Paper L 

At some later stage when the supernova shock becomes too 
weak to accelerate particles, the turbulence level in the upstream 
region goes down and the remnant can no longer hold any particles. 
At this point, all low-energy particles which remained confined un- 
til this stage escape into the interstellar medium. We assume that 
this occur when the supernova age becomes 10® yr. Then, the com- 
plete cosmic-ray escape time is taken to be 



T^sc{E) = min [tesc{E), lO'yi-] 



(4) 



Details about the nature of particle es cape fro m supernova remnants 
can be fojjnd in many literatures dPtuskin & Zirakashvili 2003, 
ICaprioli et al.l200 9. etc.). Having parameterized the cosmic-ray es- 
cape times, the coiTesponding escape radii can be calculated as. 



Resc{E) = 2.5llO tsed 



0.6 



(5) 



where uq is the initial velocity of the shock. 

Assuming that the supernova remnant is spherically symmet- 
ric, the source term in Eq. (1) is taken as 



Q{r,E,t) = 






5{t - r,sc)<5(r - Resc) 



(6) 



where Aesc{E) = 47r_Resc denotes the area of the remnant 
corresponding to the escape time of cosmic rays of kinetic en- 
ergy/nucleon E, r represents the radial variable, and q{E) — 
Aq{U) is the source spectrum with q{U) given by. 



qiU) = k{U^ + 2C/m)-(^+'^/2(f/ + m) 



(7) 



In Eq. (7), U = AE is the total kinetic energy of the particle, 
r is the source spectral index, m denotes the mass energy of the 
particle, and k is a constant. With all these ingredients, the solution 
of Eq. (1) is obtained as. 



N{rs,E,t) 



q{E) R, 



,Aesc^y'TD{t — Te, 



zexp 



R^ + r^ 

\ ^ ''esc I ' s 



X sirih 



4D(t - Tesc 

(f sJ^esc 
2D{t - Tesc 



(8) 



Eq. (8) gives the cosmic-ray spectrum at a distance Ts (measured 
from the center) from a supernova remnant injecting cosmic rays 
of different energies at different times during its evolution. Us- 
ing the properties e^ ^ 1 and sinh(a::) ~ x for very small x, 
it can be checked that at high energies where the diffusion radius 
Tdiff = \/D{t — Tesc) is much larger than {rs,Resc), Eq. (8) 
yields N{E) oc q{E)/D'"^ oc £-(r+3<5/2)^ 



where ised represents the onset of the Sedov phase, pm the maxi- 
mum rigidity of the cosmic rays and a is a positive constant. The 
maximum particle energy is assumed to scale with its charge as 
ZUm, where U,n denotes th e maximum kin etic energy for pro- 
tons. We take Um = 1 PeV JBerezhkdl 1996I) . which corresponds 



4 COSMIC-RAY BACKGROUND SPECTRUM 

4.1 Primary cosmic rays 

The background component for a primary cosmic-ray species, 
hereafter denoted by p, can be calculated from the steady state 
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diffusion-loss equation: 



V • (DpVNp) - r)Vpap5{z)Np = -Qp 



(9) 



Tfie terms on tiie left represent diffusion and losses due to inelastic 
collisions, where 77 is the averaged surface density of interstellar 
matter on the Galactic disk, Vp is the particle velocity and o-p{E) 
is the collision cross-section. We consider the diffusion region as a 
cylindrical disk bounded in the vertical directions at z — ±H and 
unbounded in the radial direction, as mentioned in section 2. For 
sources uniformly distributed in the thin Galactic disk, the source 
teiTn is represented by a delta function as Qp{r, E) = Sq(E)S{z), 
where S denotes the surface density of supernova explosion rate in 
the Galactic disk . Eq. (9) can be solved analytically as described in 
JThoudamlbOOSl) . and the solution at r = is given by 



N.i.i^)^^mf. 



2D„ 



where. 



fp = 



siah[K{H - z)] 



smh{KH) \Kcoth{KH) + ^^^ 






(10) 



X Ji{KR)dK 



In Eq. (10), Ji denotes the Bessel function of order 1, and R is the 
radial size of the supernova remnant distribution. Taking z — 0, 
Eq. (10) gives the cosmic-ray spectrum at the Earth. This is reason- 
able considering the fact that the position of the Earth (~ 8.5 kpc 
from the Galactic center) is well contained within the size of the 
source distribution which is taken to be 20 kpc, and also that the 
majority of the cosmic rays that reach the Earth are produced by 
sources within ~ 5 kpc, the size of the vertical halo boundary. It is 
important to mention that in Eq. (10), q{E) is taken to be the same 
as in the case of the local component given in the previous section, 
thus maintaining the same source spectrum between the local and 
the background components. 



4.2 Secondary cosmic rays 

Cosmic-ray secondaries are produced as spallation products from 
the interaction of heavier primaries with the interstellar matter dur- 
ing their propagation through the Galaxy. For matter distribution 
on the thin Galactic disk, the secondary production rate can be cal- 
culated as. 



Q.ir,E) 



d 



nvpNp{r,E')Siz)-—apsiE,E')dE' (II) 



dE 



where s denotes secondary species, A'p represents the primary 
number density, and daps{E,E')/dE' represents the differential 
cross-section for the production of a secondary nucleus with an en- 
ergy/nucleon E from the spallation of a primary nucleus with en- 
ergy/nucleon E' . Assuming that the energy /nucleon is conserved 
during the spallation process, the differential cross-section can be 
written as. 



_d_ 
dE 



-Ops{E,E') = cjps5{E' ~E), 



(12) 



where Ops is the total spallation cross-section of the primary to the 
secondary. Eq. (II) then reduces to 



Qs(r,E) = r)VpapsNp(r,E)5{z) 



(13) 



The propagation of secondaries follows a similar equation that 
describes their primaries as given by Eq. (9), with the source term 



replaced by Eq. (13). Their differential number density is given by 

(14) 



replaced by Eq. (1 
( lThoudaiij20"^ 



N,{z,E) = rjVpapsNp(0,E)^fs 



2D, 

where A'p(0, E) is given by Eq. (10), and 
, smh[K{H~z) 



sinh{KH) \Kcoth{KH) + ^ 



y X li{KR)dK 



J 



By taking z — 0, Eq. (14) gives the background spectrum of sec- 
ondary cosmic rays at the Earth. From Eq. (14), it can be shown 
that the ratio of the secondary to the primary densities for the back- 
ground component gives a good measure of the cosmic-ray diffu- 
sion coefficient as, 

which is a well-known result in cosmic-ray propagation studies. 



5 INTERACTION CROSS-SECTIONS AND 
INTERSTELLAR MATTER DENSITY 

5.1 Cross-sections 

In this section, we present the various interaction cross-sections and 
the interstellar matter distribution that will be used in our calcula- 
tions. For protons, the inelastic interacti on cross-section is taken 
from the simple parameterization given in lKelner et al.lCOOq : 



rjp{T) = (34.3 + 1.88L + Q.25L'^) 



1- 



Tth 
T 



mb (16) 



where T is the total energy of the cosmic-ray proton, L — 
ln(T/l TeV) and Tth = 1.22 GeV is the threshold energy for 
the production of tt" mesons. For helium a nd other heavier nu clei, 
the spallation cross-sections are taken from lLetaw et al.lll983l : 



aA{E) = ao [1 - 0.62e-^/"'2 sin (lO.gX""'^*)] 



(17) 



where A represents the mass number of the nuclei, E is the kinetic 
energy/nucleon in GeV/n, X = i5/0.001 and 



cro 



: 45^°-^ [1 + 0.016 sin (5.3 - 2.63 In A)] mb 



(18) 



Both Eqs. (16) and (17) represent a good approximation to the mea- 
sured cross-section d ata d own to sub-GeV energies. For helium, as 
suggested bv lLetaw et al.ll983L we further include a correction fac- 
tor of 0.8 in Eq. (17). 

For the secondary boron production, we consider only the 
^^C and ^^O primaries as they dominate the total boron produc- 
tion in the Galaxy. Through spallation, they produce (^^B,^"B) and 
("C,^°C) isotopes. The latter further decay into (^^B,^°B) thereby 
contributing to the production of boron. For our calculations, we 
us e the tabulated production cross-sections of these isotopes given 
in lHeinbach & Simorjll995l . 



5.2 Matter density 

Since we assume that the interstellar matter is distributed in the 
thin Galactic disk, it is more relevant for our study to determine 
the surface matter density on the disk rather than the actual number 
density. For a given distribution of atomic hydrogen in the Galaxy 
uhi {r, z), the surface density at the Galacto-centric radius r can be 
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obtained as nHi(r) — J^ nHi{r, z)dz. Similarly, for the molec- 
ular hydrogen distribution, we can calculate the surface density as 
nH2 [f) = J^ nH2 (t, z)dz. The total surface density of atomic 
hydrogen is then obtained as nnir) = nHi{r) + 2nH2 {''')■ Be- 
cause cosmic rays arriving at the Earth are mostly produced within 
a distance equivalent to the halo height H — 5 kpc, only the in- 
terstellar matter distributed within a circle of radius 5 kpc is impor- 
tant for our study. For our calculations, we use the averaged surface 
density determined for this circle. 

The distribution of atomic hy droge n is taken from 



ine distribution oi atomic hydroge n is taKen trom 
iGordon & BurtonI 1 19761 and ICox et al.l Il986l while that of the 



molecular hydrogen is taken from [Bronfman et al.l Il988l From 
these distributions, we obtain the averaged surface density of 
atomic hydrogen within our 5 kpc circle as 2.85 x 10'^" atoms 
cm^^, and that of molecular hydrogen as 1.16 x 10^" molecules 
cm~^. This gives a total averaged surface density of atomic 
hydrogen of 5.17 x 10^" atoms cm^^ which is finally used for the 
present study. In addition, it is assumed that the interstellar matter 
consists of 10% helium. 



6 CONSTRAINT ON THE CONTRIBUTION OF NEARBY 
SOURCES 

Before we proceed, we will first determine the diffusion coeffi- 
cient of cosmic rays in the Galaxy. As already mentioned in sec- 
tion 4.2, the diffusion coefficient D{p) — Do/3(p/po) can be 
determined from the secondary-to-primary ratio. For our study, 
we choose the boron-to-carbon ratio since this is the most well- 
measured and well-studied ratio. The boron-to-carbon ratio calcu- 
lated using Eq. (14) is compared with the measured data in Figure 
1 . The solid line repres ents our calcu lation, a nd the data are take n 
from HEAD (Engelman n"et alj|l99(]|) , CR N JSwordv et al 



199C), 



2011 



CREAM ( ,Ahn et ah,200^ and TRACER JObermeier e~ 
We find that choosing values of Dq — 1.55 x lO'^* cm^ s~^, 
po = 3 GV and 5 — 0.54 produces a good fit to the data. Our 
calculation takes into account the effect of solar modulation using 
the f orce field approximation with modulation parameter (j> — 400 
MV JGleeson & Axford|[l968l) . 

Both the values of Dp and 5 o btained in this work are lower 
than the values adopted in IPaper a for the case of pure diffusion 
modeB We used Do = 2.9 x 1 0^^ cm^ s~^ and S = 0.6 in Paper I . 
These values were taken from iThoudamlcOOSi which were deter- 
mined using a slightly larger value of interstellar hydrogen density, 
and based on earlier measurements before CREAM and TRACER 
data became available. In Figure 1, it can be noticed that our new 
value of D nicely agrees with the measurements up to ~ (1 — 2) 
TeV/n. 

In our model which considers the effect of nearby sources on 
the observed cosmic rays, it is important to note that nearby sources 
can produce noticeable affects mostly on the primary spectrum. The 
effect on the secondaries can be neglected. This is because for cos- 
mic rays produced by nearby sources located within ~ 1 kpc from 
the Earth, the nuclear spallation time is much longer than the propa- 
gation time to the Earth. Therefore, the primaries do not get enough 
time for spallation before they reach us. By the time they undergo 



^ For t he rest of this paper, unless otherwise stated, any compar ison with 
IPaper 11 will be always with the pure diffusion model (Model A of lPaper I) . 
Comparison for the re-acceleration model (Model B of Paper I) will not be 
presented here as the comparison will look similar to what we obtain in the 
case of the pure diffusion model. 
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Figure 1. Boron-to-carbon ratio. Data: HEAO te n ^elmann et al.lll990l) , 
CRN JSwordvetalJ 119901) . CREAM JAhn et al.l l2008h and TRACER 
iObermeier et al.l201 ll) . Thick-solid line: Our calculation. 



spallation with the interstellar matter, they have already travelled 
so far that the resulting secondary flux reaching the Earth is al- 
most negligible. A detailed calculat ion on the resulting secondary 
flux can be found in lThoudamll2008l . Thus, we assume that the sec- 
ondary cosmic rays that we observe are produced entirely by the 
background primary cosmic rays. With this assumption, the sec- 
ondary spectrum can be used to determine the contribution of the 
background cosmic rays to the observed primary spectrum. Once 
the background contribution has been fixed, we can then set a limit 
on the contribution of nearby sources since the observed spectrum 
is taken to be equal to the background plus the local component. 
Thus, the secondary spectrum can put a constraint on the contribu- 
tion of the nearby sources to the observed primary spectrum. This 
will be demonstrated in the following, taking boron and their pri- 
mary nuclei, carbon and oxygen, as an example. 

Figure 2 shows the background spectra of carbon (top) and 
oxygen (bottom) calculated using Eq. (10). The dashed lines rep- 
resent the resu lts of our calcu lations, and the data are taken 
from CREAM JAhn et al.ll2009l). ATIC jPanov et alj|2007l) . CRN 



(MuUer et al 1991'), HEAO (Engelmann et al."l990') and TRACER 
(Obermeier et al. 201 1) . The calculations assume the same source 
parameters for both elements. The source spectral index and the 
source power are chosen such that the resulting boron spectrum 
best explain the measured boron data up to ~ 200 GeV/n where 
the uncertainties in the measurements are small. The boron spec- 
trum is shown in Figure 3, where the line represents our calcula- 



tion and the measurements are taken from HEAO ( Engelmann et al. 



T99d), CRN ( ISwordvetai]|l990l) and TRACER tObermeier et al. 
2011). We find that taking the source index of F = 2.31, and the 
source power Sfc(o) ~ 4.85 x 10^* ergs Myr~^ kpc~^ produces 
a good fit to the boron data, where S represents the surface den- 
sity of the supernova explosion rate as introduced in section 4. 1 , / 
is the cosmic-ray injection efficiency in units of 10*^ ergs which 
is defined as the amount of supernova explosion energy converted 
into the primary species, and the subscript C{0) denotes carbon 
(oxygen). 

Although we did not normalize our calculations either on the 
carbon or the oxygen data, it can be seen from Figure 2 that both 
background spectra are already in very nice agreement with the 
respective measurements up to ~ 1 TeV/n. This good agreement 
between the data and the background component has already set a 
very tight constraint on the contribution of nearby sources at least 
below ~ 1 TeV/n. If we allow the maximum contribution of nearby 
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Figure 2. Carbon (top) and oxyge n (bottom) energy spectra (x i?^-^^). 



Data: CREAM JAhn et alj |2009|). ATIC jPanov et alj |2007|) . CRN 
MuUeretalJ Il99ll) HEAP )Engelmann et alj Il990h and TRACER 



Obermeier et al.l201 ll) . Thick-dashed line: Background spectrum. 
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Figure 3. Boron energy spectrum (xE"^). Data: HEAO (Engelmann et al. 
[l990), CRN iSwordyetal. 1990) and TRACER tOberm eier et al.,i2011i) . 
Thick-solid line: Background spectrum. 



sources to be ~ 10% of the total observed spectrum at 10 GeV/n 
(see next section), the carbon or oxygen injection efficiency of su- 
pernova remnants is constraint to a value of fc{0) = 0.063%, and 
the supernova surface density required to maintain the background 
is obtained to be 5 = 7.7 Myr~^ kpc^^. The latter corresponds to 
a supernova explosion rate of 0.97 per century in the Galaxy. 

It should be noted that different nearby supernova remnants 
may contribute at different energies. The constraint on the cosmic- 



ray injection efficiency that we just derived is based on the con- 
tribution of those remnants contributing at low energies taken at 
10 GeV/n. But, if we assume an equal injection efficiency for all 
the supernova remnants in the Galaxy, the same constraint can 
also apply to those nearby remnants contributing at higher ener- 
gies, thereby also putting a limit on their contribution to the ob- 
served cosmic rays. Under this constraint, the total contribution of 
the nearby supernova remnants listed in Table 1 to the observed 
cosmic rays will be calculated. In all our calculations hereafter, the 
supernova rate will be taken to be the same as given above, while 
the injection efficiencies and the spectral indices will be allowed to 
vary for different cosmic-ray species and optimized based on their 
respective data. 



7 RESULTS FOR CARBON, OXYGEN AND IRON 
NUCLEI 

In Figure 2, we can see that although the background components 
agree nicely with the data up to ~ 1 TeV/n, at higher energies 
there is a discrepancy between the data and the calculations. The 
data seem to show some excess above ~ 1 TeV/n. This excess or 
hardening in the spectrum can be explained if we include the con- 
tribution of the nearby supernova remnants as shown in Figure 4. 
In the figure, the thin-solid lines represent the total contribution of 
the nearby supernova remnants, the thick-dashed lines represent the 
background spectrum and the thick-solid lines represent the total 
nearby plus background spectrum. It can be noticed that the nearby 
remnants contribute mostly above ~ 500 GeV/n explaining the ob- 
served spectral hardening. 

In Figure 4, the thin-dashed lines show the contribution of 
the dominant supernova remnants at different energies. At ener- 
gies below ~ 300 GeV/n, the nearby contribution is dominated by 
Loopl and Monogem, while above, the main contributions come 
from Vela and G299.2-2.9. The cut-off at low energies in the case 
of Vela and G299.2-2.9 is largely due to energy-dependent cosmic- 
ray escape from the remnants with some effect of slow propagation 
at those energies, and the high energy fall-off is mainly because of 
the energy-dependent propagation effect (see section 3). 

The model parameters used in our calculation for the local 
component are discussed as follows. The value of tscd depends on 
the initial shock velocity of the supernova remnant, the initial ejecta 
mass, and the gas density of the suiTounding interstellar medium. 
Typical values fall in the range of ~ (100 — lO"') yr, and we con- 
sider tg^d — 500 yr for the present work. For the initial shock ve- 
locity, we assume a value of uo = 10^ cm/s. These values give the 
cosmic-ray escape times in the range of i^ac ~ (500 — 10"") yr, and 
Resc in the range of ~ (5 — 100) pc. The cosmic-ray escape pa- 
rameter a is kept as a free parameter, but its value is assumed to be 
same for all the cosmic-ray species. For a given cosmic-ray species, 
we take the same injection efficiency and the same source index as 
those used in the calculation for the background component. And, 
all the nearby supernova remnants considered in our study are as- 
sumed to have the same set of model parameters. 

For the results shown in Figure 4, ^c{0) ~ 2.31, fc(0) = 
0.063%, 5" = 7.7 Myr"^ kpc"^ and a = 2.4. This value of a 
gives particle escape times of t^sc = (500 — 10^ ) yr for particles of 
energies 0.5 PeV/n to 1.5 GeV/n. The present value of a is slightly 
larger than the value of 2.2 adopted in lPaper j . A larger value is re- 
quired due to the smaller value for the diffusion coefficient used in 
the present study. Physically speaking, larger a means shorter con- 
finement within the remnant, while smaller D means longer time 
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Figure 4. Carbon (top) and oxygen (bottom) energy spectra (x_B °^). 
Thick-dashed line: Background spectrum. Thin-solid line: Total nearby 
component. Thick-solid line: Total nearby plus background. Thin-dashed 
lines: Dominant nearby supernova remnants. 
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Figur e 5. Iron energy spectrum {xE^'^) . Data: CREAM (Ahn et alj 



120091). ATIC JPanovetalJ |2007|) . CR N JMuUer et alj Il99ll) HEAO 
I'Enge lmann et alj Il990l) ~and TRACER iObermeier et al.l l201ll) . Thick- 
dashed line: Background spectrum. Thin-solid line: Total nearby compo- 
nent. Thick-solid line: Total nearby plus background. Thin-dashed lines: 
Dominant nearby supernova remnants. 



Table 2. Source spectral indices F and injection efficiencies / for the vari- 
ous cosmic-ray nuclei considered in our study. 



Nuclei 


r 


/( 


X 10"*^ ergs) 


Proton 


2.27 




17.5 


Helium 


2.21 




1.75 


Carbon 


2.31 




0.063 


Oxygen 


2.31 




0.063 


Iron 


2.31 




0.011 



for cosmic rays to reach the Earth. Therefore, in order to explain 
the spectral hardening above ~ 500 GeV/n, which in our model is 
due to the contribution of the nearby sources, a larger a is required 
to compensate the effect of a smaller D. This will become more 
clear in the next section, when we compare our present r esults for 
protons and helium nuclei with those obtained in lPaper B . 

For iron, the result is shown in Figure 5. We take the same 
source index of Tpe = 2.31 as in the case of carbon and oxygen, 
and an injection efficiency of /fe = 0.011%. All other model pa- 
rameters remain the same as in Figure 4. Also, in the case of iron, 
we can see that the nearby supernova remnants produce a spectral 
hardening above ~ 500 GeV/n. Future sensitive measurements at 
high energies can provide a crucial check of our prediction. 



8 RESULTS FOR PROTON AND HELIUM NUCLEI 

In Figure 6, we show our results for protons (top) and helium nu- 
clei (bottom). The thick-dashed and the thin-solid lines represent 
the background and the total nearby contributions respectively, and 
the thick-solid line shows the total sum of the background and the 
nearby components. The contributions from the nearby dominant 
sources are also shown by the thin-dashed lin es. The m easurements 
data are from CR EAM JYoon et al.l 1201 ll) . ATIC0 franov et al.' 
l2007h , PAMELA^ JAdriani et al.ll201lh , and AMS ( Alcaraz et al. 



Data taken from the compilation bv lStrong & Moskalenkoll2009l 



I2OOCI : lAguilar et al.ll2002l) . The results look very similar to those 
obtained for the heavier nuclei in section 7. Except for the source 
index and the injection efficiency, all other model parameters are 
taken to be the same as for the heavier nuclei. For protons, we find 
that taking Tp = 2.27 and fp = 17.5% produces a good fit to 
the data, while for helium the best fit parameters are found to be 
Vhb = 2.21 and fne ~ 1.75%. For comparison, these values are 
listed in Table 2 along with those obtained for carbon, oxygen and 
iron. 

In gene ral, our present results are quite similar to the results 
presented in lPaperli For instance, the nearby supernova remnants 
show dominant contribution at energies > (0.5 — 1) TeV/n, thereby 
explaining the observed spectral hardening. Moreover, the helium 
spectrum show hardening at lower energy/nucleon with respect to 
the proton spectrum which, in our model, is attributed mainly to 
the early escape times of helium nuclei from the supernova rem- 
nants relative to the protons. It might be recalled from lPaper j. and 
also discussed in section 3 of this paper that in our model, such 
a spectral hardening at lower energies/nucleon is expected for all 
heavier primaries {A > 1) whose escape times are shorter than 
the time for protons by a factor of (A/Z)"^'" at the same en- 
ergy/nucleon. Also, both our present and the previous studies show 
that the main contribution at high energies comes mainly from the 
Vela and G299.2-2.9 remnants, and that the spectral hardening does 
not continue up to very high energies. Available measurements also 
seem to support these results. 

However, some basic differences can be noticed. First, the 
background spectrum is steeper in the present case. At high en- 
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Figure 6. Proton ( x_E^'^5^ i-^p^ ^^^j helium ( x_E^-^, bottom) e nergy spec- 
tra. Data: CREAM JYoon et alj201lh. ATIC jPanov et alj2007h. PAMELA 



JAdriani et ai]l201 11) . and AMS JAlcaraz et al.ll2OO0l: lAguilar et al.ll2002l) . 
Thick-dashed line: Background spectrum. Thin-solid line: Total nearby 
component. Thick-solid line: Total nearby plus background. Thin-dashed 
lines: Dominant nearby supernova remnants. 



ergies, the proton background follows a spectral index of 2.82 and 
the helium background has an index of 2.75. To be compared, the 
b ackgroun d indices were obtained as 2.75 and 2.68 respectively 
in lPaper j . It may be recalled that in jPaper B, the backgrounds were 
obtained by fitting the measured spectra between (20 — 200) GeV/n 
and their consistency with the low energy data below ~ 20 GeV/n 
was not checked. From Figure 6, it can be noticed that the back- 
ground adopted in the present study agrees nicely (in fact, even 
better when the small local component has been added) with the 
data down to 1 GeV/n. 

Second, in IPaper 11 the source index for the nearby sources 
were obtained as (Fb — 5), where F;, is the background index and 5 
is the diffusion index. S o, for the proton background of Fj, = 2.75 
and 5 = 0.6 adopted in Paper J. the source index was found to be 
2.15. This is flatter than the proton source index of 2.27 adopted in 
the present study. Adopting a steeper source index suppresses the 
contribution from nearby sources at high energies. For the same 
amount of total energy injected into protons, a source with an in- 
dex of 2.27 produces '--^ 1.8 times less number of source particles at 
~ 10 TeV than a source with an index of 2. 15. A similar difference 
is also expected in the case of helium. On the other hand, taking a 
smaller value of D enhances the flux from a nearby source. This 
is clear from the discussion on Eq. (8), given in section 3 which 
showed that at very high energies, the particle spectrum depends 
on D as TV oc D~^'^ . After detailed investigation, we find that the 



increase in the particle flux in the present study due to smaller D 
is almost equal to the decrease in the flux due to the steeper source 
index. Because of these two almost equal and opposite effects, we 
can still explain the spectral hardeni ng of he lium with an injection 
efficiency very close to that used in lPaperj . For protons, we need 
an injection efficiency of 17.5% in the present s tudy, which is ap- 
proximately twice the value obtained in lPaperll This difference is 
because of the lower proton background in the present stud y at en- 
ergies above ~ 1 TeV as compared to the background in Paper B 
For helium, the background does not differ too much at TeV ener- 
gies in the two studies. 

Another difference results from the difference in the cosmic- 
ray escape parameter and the diffusion index. A smaller a produces 
a sharper low-energy cut-off and a larger 5 leads to a stee per fall- 
off in the high-energy spectrum from a nearby source. In IPaper ll 
where we took a = 2.2 and 5 = 0.6, this led to sharper peaks 
in the individual contributions of the nearby remnants, thereby re- 
sulted into stronger structures in the resultant total spectrum. In the 
present study, the slightly larger value of a = 2.4 and the smaller 
value of (5 = 0.54 produce broader peaks in the individual contri- 
butions leading to weaker structures in the overall total spectrum. 



9 DISCUSSION 

We have shown that the spectral hardening of cosmic rays at TeV 
energies recently observed by the ATIC, CREAM, and PAMELA 
experiments can be due to nearby supernova remnants. Consider- 
ing that cosmic rays escape from supernova remnants in an energy- 
dependent manner, we also show that heavier elements should pro- 
duce spectral hardening at relatively lower energies/nucleon with 
respect to the protons, and that the hardening might not continue 
up to very high energies. These results are also seemed to agree 
with the measured data. 

In general, the res ults obta ined in this paper agree very well 
with those presented in IPaper B . Our present study involv es a de- 
tailed calculation of the background cosmic rays unlike in lPaperB, 
and also follow a consistent treatment of the cosmic-ray source 
spectrum for the background and the nearby sources. Our results 
are found to be consistent with the observed data over a wide range 
in energy from 1 GeV/n to ~ 10^ GeV/n for a reasonable set of 
model parameters. Our calculation requires a supernova explosion 
rate of ~ 1 per century in the Galaxy, and cosmic-ray injection ef- 
ficiencies of fp = 17.5% for protons, /hs = 1.75% for helium 
nuclei which is exactly 10% of the proton value, fc(0) = 0.063% 
for carbon and oxygen, and fpe = 0.011% for iron. The re- 
quired source index for protons is Fp — 2.27 and for helium nu- 
clei, Tue = 2.21. For carbon, oxygen and iron, we determined 
the same source index of 2.31. The required source indices of 
~ (2.2 — 2.3) in the present study are slightly steeper than the 
value of F = 2.0 — 2.2 predicted by DSA theory. Actually, even 
a larger value of F ~ (2.4 — 2.5) is favored by the high level 
of cosmic-ray isotropy observed between around 1 and 100 TeV 
which in t urn suggests a sma ller value of the diffusion index of 5 ~ 
0.3 to 0.4 ( jPtuskin et alj2006.1 . This discrepancy between observa- 
tion and theory is still not clearly understood. 

Our model predictions are expected to be different in many re- 
spects from those of other existing models. Models which are based 
either on the hardening in the source spectrum or changes in the 
diffusion properties of cosmic rays at high energies will produce 
a spectrum that remains hard up to ver y high energies (lOhira et al.l 
l201lllYuan et al.l20IlllAve et al.l20o"^ . But, although not very sig- 
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Figure 7. Proton energy spectrum under different models. Thick-dashed 
line: Model I. Thin-solid line: Model II. Double-dashed line: Total back- 
ground plus nearby spectrum in our model. Thin-dashed line: Our back- 
ground spectrum. Data as given in Figure 6 (top). 
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Figure 8. Boron-to-carbon ratio under different models. Thick-dashed line: 
Model I. Thin-solid line: Model II. Double-dashed line: Our model. Data as 
given in Figure 1 . 



nificant, the CREAM data seems to indicate that the spectral hard- 
ening does not continue beyond few tens of TeV/n, which in general 
is in good agreement with our prediction. In Figure 7, we compare 
our result for proton spectrum (double-dashed line) with the pre- 
dictions of other models: Model I (thick-dashed line) and Model II 
(thin-solid line). Model I represents models with a hardened source 
spectrum above a certain energy, but assumes the same diffusion 
coefficient as in our model. Model II represents models which in- 
corporate a break or hardening in the diffusion coefficient. The 
source spectrum in Model II is taken to be the same as in our model. 
Also shown in Figure 7 for reference is the background spectrum 
obtained in our model. To reproduce the same measured spectrum, 
Models I and II are chosen to have their respective breaks in the 
source spectrum or in the diffusion coefficient at the same energy 
Eb = 850 GeV/n. It can be noticed that our result shows significant 
difference mostly at energies above '--^ 0.1 PeV. Thus, if the cosmic- 
ray spectrum exhibit an exponential cut-off below '^ 0.1 PeV, our 
result will not be significantly different from the others. However, 
detailed studies on the origin of the cosmic-ray knee suggests a 
cut-off at energies around 3 — 5 PeV (H orandel 20031) . 

Even more different between the different models will be the 
secondary-to-primary ratios and the secondary spectra. In the stan- 
dard model of cosmic-ray propagation, the secondary-to-primary 
ratio is independent of the source parameters and gives a good mea- 
sure of the cosmic-ray diffusion coefficient in the Galaxy. However, 
in our model, which considers the effect of the nearby sources, 
the ratio may deviate from the standard result. This is because the 
nearby sources can affect only the primary spectrum, and their ef- 
fect on the secondaries is negligible. Thus, we expect to see a steep- 
ening in the ratio in the energy region where the nearby contribu- 
tion on the primary spectrum is significant. This is shown in Figure 
8 for the boron-to-carbon ratio. Notice the significant steepening 
in our model above ~ 500 GeV/n with respect to Model I even 
though both the models assume the same diffusion coefficient. The 
result for Model II is even more flatter at high energies, reflecting 
the harder value of diffusion coefficient above Eb — 850 GeV/n 
assumed in the model. 

For an equilibrium primary spectrum Np oa E^"' , and a dif- 
fusion coefficient D oc E , the secondary spectrum in the Galaxy 
follows Ns oc i5~ '^^ ' . This shows that for a fixed diffusion coef- 
ficient, the shape of the secondary spectrum is determined by the 
shape of the primary spectrum. Therefore, models that consider 
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Figure 9. Boron energy spectrum under different models. Thick-dashed 
line: Model I. Thin-solid line: Model II. Double-dashed line: Our model. 
Data as given in Figure 3. 



the same diffusion coefficient, but different Np will produce dif- 
ferent A'^s. Since the background primary spectrum above Eb in 
our model is steeper than in Model I (see Figure 7, thin dashed 
and thick-dashed lines), the secondary spectrum is expected to be 
steeper above Eb in our model. This is shown in Figure 9 for the 
boron spectrum, where the thick-dashed line represents Model I, 
and the thin-dashed line represents our model. The difference is 
expected to be even more significant if we compare with Model II, 
which assumes a harder diffusion index above Eb. The thin-solid 
line in Figure 9 represents Model II. Similar differences are also 
expected in other types of seconda ry nuclear species such as sub- 
iron, and anti-protons as shown in Iviadimirov et al .1120121 which 
also discussed the effects of different models on various observed 
properties of cosmic rays including the secondary-to-primary ratios 
and the diffuse gamma-ray emissions. 

The Galactic diffuse gamma-ray spectrum in our model is also 
expected to be different from those calculated using other models. 
If the diffuse emission is dominated by gamma-rays produced from 
the decay of tt" mesons, then their spectrum at high energies would 
largely follow that of the primary protons. As the background spec- 
trum in our model is steeper above Eb than the spectrum adopted 
in other models, our diffuse gamma-ray spectrum will be steeper 
above Eb. This is shown in Figure 10 where the thin-dashed line 
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Figure 10. Galactic diffuse gamma-ray spectrum under different mod- 
els, nick-dashed line: Model I. Thin-solid line: Model II. Double-dashed 
line: Our model. The data represents the gamma-ray intensity measured by 
FERMI for Galactic latitudes \b\ ^ 10° l Abdo et al. 2010) . 



represents our model, and the thick dashed and thin-sohd Unes rep- 
resent Models I and II respectively. The data represents the gamma- 
ray intensity for Ga lactic latitudes |6 | Js 10° measured by the 
FERMI experiment JAbdo et al]|201(3h . All results plotted in Fig- 
ure 10 are normalized to the data at 1.2 GeV. The calculations use 
the different proton spectra shown in Figure 7, and the gamma-ray 
production cross-section is taken from Kelner et al. 2006. In Figure 
10, we can see that Models I and II give almost the same result, and 
produce a harder gamma- ray spectrum with respect to our model 
above ~ 100 GeV. This difference can be checked by future mea- 
surements at high energies, and can distinguish our model from 
others. The data in Figure 10 show some excess above the model 
predictions between '--^ (10 — 100) GeV. Although it is not the aim 
of this paper to perform a detailed modeling of the FERMI data, it 
can be mentioned that the excess is most likely due to additional 
contributions from other processes such as bremsstrahlung, inverse 
compton and unresolved point sources, which are neglected in our 
calculations. Detailed calculations involving all the possible con- 
tributions have shown that the diffuse gamma-rays up to --^ 100 
GeV measured by FERMI from different regions in the Galaxy can 
be explained using a single powe r-law cosmic-ray spectr um with- 
out any break above a few GeV JAckermann et al.ll20I2l see also 
IVladimirov et al.l2012h . 



10 CONCLUSION 

We have presented a detailed and improved version of our previ- 
ous work presented in IPaper 11 where we showed that the spectral 
hardening of cosmic rays observed at TeV energies can be a local 
effect due to nearby supernova remnants. Unlike in Paper 11 where 
the cosmic-ray background was obtained by merely fitting the low- 
energy data, the present work involves a detailed calculation of both 
the background and the local components considering consistent 
cosmic-ray source parameters between the two components. 

In addition to th e results for protons and helium nuclei (which 
were also studied in IPaper 11) ■ we have also presented results for 
heavier cosmic-ray species, such as boron, carbon, oxygen, and 
iron nuclei. Unlike in other existing models, we have shown that 
heavier nuclei should exhibit hardening at lower energy/nucleon 
as compared to protons, and also that the hardening may not con- 
tinue up to very high energies for all the species. Although not 



very significant, the available data seem to suggest our findings. 
Moreover, we have also shown that our results on the secondary- 
to-primary ratios, the secondary spectra, and the Galactic diffuse 
gamma-ray spectrum at high energies are expected to be different 
from the predictions of other models. Future sensitive high-energy 
measurements on these quantities can differentiate our model from 
others. 
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